Local moment formation and Kondo screening in impurity trimers 



Andrew K. Mitchell, Thomas F. Jarrold, Martin R. Galpin, and David E. Logan 

Department of Chemistry, Physical and Theoretical Chemistry, 
Oxford University, South Parks Road, Oxford 0X1 3QZ, United Kingdom 
(Dated: March 5, 2013) 

We study theoretically a triangular cluster of three magnetic impurities, hybridizing locally with 
conduction electrons of a metallic host. Such a cluster is the simplest to exhibit frustration - 
an important generic feature of many complex molecular systems in which different interactions 
compete. Here, low-energy doublet states of the trimer are favored by effective exchange interactions 
produced by strong electronic repulsion in localized impurity orbit als. Parity symmetry protects a 
level- crossing of such states on tuning microscopic parameters, while an avoided crossing arises in 
the general distorted case. On coupling to a metallic host, the behavior is shown to be immensely 
rich, since collective quantum many-body effects now also compete. In particular, impurity degrees 
of freedom are totally screened at low temperatures in a Kondo-screened Fermi liquid phase, while 
degenerate ground states persist in a local moment phase. Local frustration drives the quantum 
phase transition between the two, which may be first order or of Kosterlitz-Thouless type, depending 
on symmetries. Unusual mechanisms for local moment formation and Kondo screening are found, 
due to the orbital structure of the impurity trimer. Our results are of relevance for triple quantum 
dot devices. The problem is studied by a combination of analytical arguments and the numerical 
renormalization group. 

PACS numbers: 



I. INTRODUCTION AND MOTIVATION 

The interplay of orbital and spin degrees of freedom 
with electronic interactions can produce a diverse range 
of chemical and physical behavior. At the few-electron 
level of a single molecule, understanding the resulting 
complexity is a traditional problem in theoretical chem- 
istry — and a challenge because strong correlations pre- 
clude an independent particle picture. At the many- 
electron level of clean bulk metals, by contrast, interac- 
tions are often rather unimportant, and the system ade- 
quately described by an essentially independent particle 
description. 

Bridging between these limits is the fascinating class 
of 'quantum impurity problems', 1 in which an interacting 
and, in effect, small molecular subsystem is coupled lo- 
cally to the continuum of conduction electrons in a metal. 
The most basic example, a classic paradigm in condensed 
matter science, 1 3 is a single magnetic impurity embed- 
ded in a metallic host. Local Coulomb repulsion favors 
single-occupancy of the active impurity orbital, with a 
local moment thus forming at high temperatures due to 
the free spin degree of freedom. But the Kondo effect 
plays a key role at low temperatures/energies: the im- 
purity moment is screened dynamically by conduction 
electrons, which together form a many-body spin-singlet 
state (the 'Kondo singlet'). The screening process itself 
can be understood in terms of a renormalization group 
flow, corresponding to the crossover from the local mo- 
ment fixed point applicable at high temperatures, to the 
strong coupling fixed point describing the Kondo sin- 
glet ground state. 1 ' 4 ' 5 Importantly, physical properties 
exhibit universality in terms of the crossover tempera- 
ture scale, Tk- 

The Kondo effect associated with screening of a sin- 



gle magnetic impurity has been observed experimen- 
tally in many metals. 1 In semi-metallic systems by con- 
trast, where the host density of states vanishes at the 
Fermi level, it is well known that Kondo physics can 
be suppressed, and that such systems also support de- 
generate local moment ground states. 6 10 Local quan- 
tum phase transitions between Kondo screened phases 
and local moment phases have thus attracted much in- 
terest and have been studied in detail; being sought for 
example in softgapped systems such as graphene, 11 ' 12 d- 
wave superconductors 13 and surfaces of 3d topological 
insulators. 14 ' 15 

The situation is unsurprisingly richer when several 
impurities are present. Already in the case of two 
impurities, the resulting behavior can be markedly 
different. 16 26 If the impurities are spatially well sep- 
arated, each is essentially screened independently by 
the Kondo effect. But an effective exchange interac- 
tion ('RKKY') can arise and dominate when two impu- 
rities are brought closer, with a resulting crossover to lo- 
cal interimpurity spin-singlet formation 16 ' 22 (indeed for 
strictly independent screening channels, the crossover is 
sharpened to a true transition 17 ' 19 21 ). 

Multi-impurity systems, and even small molecules de- 
posited on metallic surfaces, exhibit more complex phe- 
nomena due to competition between local interactions 
and the Kondo effect. 27 37 In this paper we consider a 
cluster of three impurities, which exhibits the frustrat- 
ing effects of competing ground states. Frustration and 
degenerate ground states at level crossings are of course 
an important general feature of many complex molecular 
systems familiar in chemistry. Here we examine the in- 
terplay between frustration and the Kondo effect in one 
of its simplest realizations, to obtain a detailed under- 
standing of the underlying physical behavior. 
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protected level crossing. 29 In the more general distorted 
case by contrast, the transition occurs at the critical 
end point of a line of Kondo screened states. An effec- 
tive model is derived to describe this Kosterlitz-Thouless 
transition, 49 which captures the characteristic vanish- 
ing of the low-energy scale, T#, as the transition is ap- 
proached from the Kondo phase. Analytic arguments are 
supplemented and confirmed by exact numerics obtained 
using the numerical renormalization group. 4 



FIG. 1: Schematic illustration of the TQD setup. Transport 
measurements between source (s) and drain (d) leads indi- 
cated by the arrow. 



Real three-impurity systems, such as Cr clus- 
ters on clean Au(lll) surfaces or Co clusters on 
Cu(100), have been studied both experimentally 27,31 and 
theoretically. 28,32 ' 38 40 But since the impurities couple to 
different surface sites of the host, a realistic theoretical 
model maps irreducibly to a three-channel problem, 28 
whose properties depend quite sensitively on details of 
the particular experimental realization. More general 
aspects of the three impurity problem can however be 
accessed in triple quantum dot devices, 41 44 where the 
metallic leads to which the dots are tunnel-coupled pro- 
vide the conduction electrons. The Kondo effect has been 
observed in single semiconductor quantum dots coupled 
to metallic leads, 45 which are often referred to as 'arti- 
ficial atoms'. 46 By the same token, coupled dot devices 
behave as highly tunable 'artificial molecules'. 42,47 ' 48 The 
primary experimental probe in all cases is of course 
the conductance, obtained by driving a steady current 
through the lead-coupled quantum dot system. 

The system studied here is illustrated in Fig. 1, and 
consists of a triangular cluster of three quantum dots 
(impurities), tunnel-coupled apically to a single channel 
of conduction electrons (and with the metallic lead split 
into source and drain to allow measurement of conduc- 
tance). We find both a Kondo screened phase and a free 
local moment phase, depending on model parameters (see 
Fig. 2), which themselves could be tuned in a real device 
by application of gate voltages. 41 44 The Kondo phase is 
characterized by a strong coupling Fermi liquid ground 
state, in which all dot degrees of freedom are screened. 
Zero-bias conductance through the device in this case 
is maximally enhanced at low temperatures due to the 
Kondo effect. Three distinct screening mechanisms are 
uncovered, depending on the relative strength of inter- 
dot tunnel-couplings, which have characteristic univer- 
sal signatures in physical quantities. By contrast, the 
Kondo effect is totally suppressed in the local moment 
phase, which has a doubly-degenerate ground state and 
extremely low conductance. 

The quantum phase transition separating the Kondo 
and local moment phases is studied in some detail. In 
the parity-symmetric case, the transition is a parity- 



II. MODELS AND METHODS 

We consider the triple quantum dot (TQD) device il- 
lustrated in Fig. 1. It consists of three equivalent and 
locally correlated single-level ('Andersonian') sites, with 
level energy e and on-site Coulomb repulsion U. Dots i 
and j are tunnel-coupled by a matrix element to form 
a triangular arrangement. Dot 1 is also coupled to source 
and drain leads (and we consider the zero-bias case where 
the system is in equilibrium). The Hamiltonian is decom- 
posed as H = i^ieads + #tqd + #h y b where, in standard 
notation, 
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with hi a = d\ a d io . the number operator for spin-cr =f / j. 
electrons on dot site i = 1,2,3. Htqb describes the iso- 
lated TQD, by itself a small quasi- molecular entity that 
can accommodate up to 6 electrons. In the absence of 
electron interactions (f7 = 0) it is in fact merely a 3-site 
Hiickel ring, although interactions are essential, since U 
is inversely proportional to the dot capacitance (and no 
quantum dot has infinite capacitance!). i?i ea ds by con- 
trast describes the leads, an essentially non-interacting 
but macroscopic metal. The hybridization term i^hyb, in 
connecting the two subsystems via tunnel-coupling be- 
tween dot 1 and the leads, ensures that the lead-coupled 
TQD system is an interacting, many-body problem con- 
taining macroscopic numbers of electrons. That is essen- 
tial to the basic physics (despite some naive attempts to 
avoid it in the literature). 

For equivalent leads a = s, d, we define the conven- 
tional 'local' orbital to which dot 1 couples as, 
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(2) 



With this, the full TQD problem can then be mapped 
to an effective Id problem by tridiagonalizing the con- 
duction electrons, 1 ' 4 starting from the fo a 'zero' orbital. 
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One thereby obtains, 



-Hleads =^ ^ [en/ncr/ncr + h n {flaf(n+l)a + H.C.)] 

(3a) 



(3b) 



with all diagonal one-electron energies e n = if the local 
conduction electron density of states (DOS), p£\uo), is 
particle-hole symmetric; with p^\oo) given by 



^ ) H = E(^)"^- e ^) = -'imG^M (4) 



7T 



where Gq*J(uj) = ((/ 0(T ; /cL))^ * s ^ ne ^ ree Green function 
for the local bath site, fo a (correlation functions of the 
type ((A^B)}^ are simply the Fourier transform of the 
usual retarded functions —iO{t\ — t2){{A(t\) 1 Bfa)})). 

The conduction electron and hybridization terms of the 
Hamiltonian enter in such quantum impurity problems 
through a single 'hybridization function', 1 given generally 

by 



r a (w) = -v 2 im GtV) = r?M - irj( w ) 



(5) 



with r^.(cj) = ttV 2 p^\uo) , and where the real part, 
r^(cj), follows by Hilbert transformation of the imag- 
inary part, r^(cj). As such, the free DOS and the 
tunnel-coupling V specify completely the effect of the 
leads on the TQD. For simplicity, we take a wide flat 
DOS, which is the conventional and generic case 1 rele- 
vant to most metallic hosts: p^\uo) = ^0(D — 
such that T^(cj) = T = irpV 2 is defined inside the band, 
-D < uo < D, and p = 1/(22}). 



A. Numerical Renormalization Group 

The Numerical Renormalization Group 4 (NRG) is a 
non-perturbative technique for treating quantum impu- 
rity problems, such as the lead-coupled TQD system (for 
a recent review, see Ref. 50). Numerically exact thermo- 
dynamic and dynamic quantities can be calculated with 
NRG on essentially all relevant energy scales — from the 
conduction electron bandwidth D (typically a few eV in 
bulk metals), down to the Kondo scale Tk (which is an 
exponentially-small 1 fraction of D). 

The first step 4 is to divide the conduction electron DOS 
into intervals whose width decreases exponentially (the 
points separating intervals are uo n = ±DA _n , with A > 1 
and n = 0,1,2...). The spectrum is then discretized by 
replacing the DOS in each interval by a single pole of the 
same total weight. The 'Wilson chain' 4 ' 5 is then defined 
by writing the conduction electron Hamiltonian as a Id 
tight-binding chain of form Eq. (3a), and parameters /i n , 



e n chosen so that the free DOS at the TQD site cor- 
responds to the discretized spectrum. This discretized 
model is then diagonalized iteratively, with high-lying 
states discarded at each step (the truncation being jus- 
tified because the strength of the hoppings h n ~ A -n / 2 
decrease exponentially down the chain 4 ' 5 ). 

One thus examines the behavior of the system on pro- 
gressively lower energy scales, and thermodynamics can 
be built up as a function of temperature. 4 ' 5 ' 50 Dynam- 
ical quantities, such as spectral functions, can also be 
obtained from the full density matrix, 51 ' 52 calculated it- 
eratively in the Anders- Schiller basis. 53 For the calcu- 
lations presented in this paper, we use a discretization 
parameter A = 3, and retain around 3000 states at each 
iteration. 



B. Physical quantities 

We consider below the contribution to thermodynam- 
ics arising from the lead-coupled TQD system, specifi- 
cally the so-called excess quantities 4 ' 5 (0)i mp = (Cl) — 
(^)o, with (O)o denoting the thermal average in the 
absence of the TQD itself. We focus on the en- 
tropy Si m p(T), uniform spin susceptibility Ximp(^) = 
((S z ) 2 )imp/T, and 'excess' charge n imp = (Q) imp (here 
S z and Q refer to the spin and charge of the entire sys- 
tem). These quantities show characteristic signatures of 
the underlying fixed points, reached under renormaliza- 
tion on progressive reduction of the temperature/energy 
scale. The renormalization group flow between such fixed 
points shows up as crossover behavior in thermodynam- 
ics, which thus allows determination of emergent energy 
scales in the problem, such as the Kondo temperature, 

We also consider dynamical properties — in partic- 
ular dot spectral functions (or, loosely, local densities 
of states), given by D ia (uo) = — ^Im Gu^{uo) in terms 
of the local Green function for that dot, Gu^{uo) = 
((d ia ; <$ i(7 ))u)' I n the zero-bias (equilibrium) limit of inter- 
est, it is naturally the local spectrum of dot 1 (itself di- 
rectly tunnel-coupled to the leads) which determines the 
differential conductance G C (T) between source and drain 
leads, mediated via the TQD. This is given exactly 54 by 



Gc(T) 
Go 



dm 

duo 



7rTD 1(J (uo,T) duo 



(6) 



where f(uo) = [1 + exp(cj/T)] _1 is the Fermi function 
(and uo = is the Fermi level). Here, Go = 2e 2 h~ 1 Go, 
where the dimensionless quantity Go = 4F s Td/(T s +T^) 2 
embodies simply the relative coupling strength of dot 1 to 
source and drain leads; such that Go = 2e 2 h~ 1 is maxi- 
mal in the symmetric case T s = T^. At zero temperature, 
the conductance is thus controlled by the behavior of the 
spectrum at the Fermi level, 



Ge(0) 

Go 



= irTD la (w = 0, T = 0) = sin 2 (S) , (7) 
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FIG. 2: Schematic phase diagram for the TQD in the parity-symmetric case, with e = —U/2. For a given U (^> T), both 
local moment (LM) and Kondo strong coupling (SC) phases can be accessed on tuning the interdot tunnel-couplings t and t' . 
Regimes and crossovers discussed in the text. The QPT between LM and SC phases indicated by arrows (a)-(c) is considered 
explicitly in the following sections. 



where the final identity follows from the definition of the 
scattering phase shift 5 = arg[Gi cr (co' = 0,T = 0)], to- 
gether with the fact that all electron scattering vanishes 
at the Fermi level. 



III. RESULTS: SYMMETRIC CASE 

We consider first the full TQD model, Eq. (1), in the 
parity-symmetric case ti2 = tys = t and £23 = t' (see 
Fig. 1). In this limit the Hamiltonian is invariant to 
swapping the dot labels 2 <H> 3, meaning formally that 
[H, P23] = for the permutation operator P23- All states 
can thus be classified according to parity P23 = ±1 (since 
P23 = !)• As shown in Ref. 29 this allows for the pos- 
sibility of a (zero temperature) quantum phase transi- 
tion (QPT) between two parity-distinct phases, which 
are of Kondo screened and local moment types. That 
scenario is explored in detail below, focusing primarily 
on the strongly correlated case U ^> T. Since the most 
interesting physical behavior arises when the TQD is in 
essence triply occupied, we consider the representative 
case e = —U/2. The largest energy scale of the problem, 
the bandwidth D, is taken as D = 100T. 



A. Phase diagram 

Before discussing specific results for the symmetric 
TQD, in Fig. 2 we present a schematic phase diagram 
which highlights the relevant regimes, crossover scales 
and phases. The parity-symmetric TQD model supports 



two phases: a Kondo screened phase with strong coupling 
(SC) ground state, and a local moment (LM) phase with 
a doubly-degenerate ground state. 29 

As shown below, however, the effective spins which re- 
main free in the LM phase depend on the underlying 
parameters in the (£, t') -plane. Likewise, Kondo screen- 
ing proceeds by different mechanisms in the SC phase 
depending on relative magnitudes of the interdot tunnel- 
couplings. When the interdot tunnel-couplings are large 
compared to the interaction, t ^> U (regime (a) of Fig. 2), 
the relevant TQD states are simple molecular orbitals 
(MOs); with the lowest doublet state remaining free down 
to T = when t' < U. By contrast, for t,t' <C U each 
dot is essentially singly-occupied. In this case, dot 1 (con- 
nected directly to the leads) can undergo the Kondo ef- 
fect on the scale of Tk — provided the interdot tunnel- 
couplings are weak (t 2 /U <C Tk)- The strength of the 
effective coupling between the remaining dots 2 and 3 
then determines the ultimate ground state of the system, 
with a level-crossing QPT occurring between LM and SC 
states as t' is increased through t' ~ t 2 / \JUTk (regime 
(c) of Fig. 2). If by contrast the interdot couplings are 
stronger, t 2 /U ^> Tk, the lowest-energy TQD doublet 
state has weight on all three dots. This results - regime 
(b) of Fig. 2 - in an effective antiferromagnetic coupling 
to the leads if t' > t, in which case the Kondo effect is 
operative and the SC phase obtains; or, if t' < t, in an 
effective ferromagnetic coupling to the leads, such that 
the TQD doublet remains free down to T = in this LM 
phase. 

In the following, we discuss in some detail the behavior 
in these three distinct regimes. 
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B. Molecular orbital regime 

We examine first the simplest case, in which the in- 
terdot tunnel-couplings are strong, t ^> U. Here, the 
isolated TQD states are essentially non-interacting MOs. 
To obtain a handle on the problem in this regime, we 
thus derive a low-energy effective model for the non- 
interacting TQD (U = 0), then incorporate the effect 
of interactions within lowest-order perturbation theory 
mU. 

When U = 0, a simple canonical transformation of 
the dot operators {d 1(J , d 2a , d 3a } ->> {d ua , d i(J , d oa } brings 
Htqb in Eq. 1 to diagonal form, 

H tqb = J2( Euhu ° + E ° h °" + ( 8 ) 

where h aa = d^ aa d aa (a = u,l,o). The single-particle 
levels (MO energies) are given by E Q = e — t' and E u /i = 
e + \t' ± \\J%t 2 + t f2 , and for t' < t have the relative 
ordering E\ < E Q < < E u . The difference in energy 
between the 4-electron state and the 3-electron state is 
thus E/\ = E Q = e — t f < 0. The lowest many-particle 
TQD state is thus always a 4-electron state when U = 0. 

Intuition suggests however that the Coulomb repul- 
sion will favor a 3-electron state. Treating perturba- 
tively the interaction part of the Hamiltonian, Hi = 
U^2 i hi^hi^ we thus calculate the correction to the iso- 
lated TQD states to first-order in U. The resultant en- 
ergy difference between the 4- and 3-electron states is 
then E' A = e-t' + \U{1 + 7 2 ) + 0(U 2 ), with 2 7 2 = 
1 - [1 + 8(t/0 2 ]" 1/2 - Hence for e = -U/2 as considered 
explicitly, E A = U/A — t', indicating that the 3-electron 
TQD state is indeed favored when t' <U/4. 

Of course, the interesting behavior arises on coupling 
the TQD to the leads. In the relevant case of T <C U 
(<C t), one can obtain an effective model valid at low tem- 
peratures/energies <C E A , by projecting the full Hamil- 
tonian onto the 4-electron TQD state (for E A < 0) or 
the 3-electron state (for E A > 0). This is achieved by a 
Schrieffer- Wolff transformation (SWT) 55 to second order 
in i^hyb, eliminating perturbatively excitations to higher- 
lying TQD states. 

In the effective 4-electron sector, the spin-singlet TQD 
state essentially decouples from the leads. This 'frozen 
impurity' fixed point is continuously connected to the 
Kondo screened ground states 1 ' 5 (discussed further be- 
low). More interestingly, in the 3-electron sector of the 
TQD, one obtains via a SWT the effective model 

H 6ff = Pleads + 4S-S(0), (9) 

with S a spin- ^ operator for the lowest (doublet) 3- 
electron TQD state, and s(0) the spin density of con- 
duction electrons at the TQD. Eq. 9 is a model of Kondo 
form, 1 with its underlying physics well known to depend 
crucially on the sign of the exchange coupling, j£. Im- 
portantly, we find here j£ < 0, wherein the ferromag- 
netic Kondo effect arises. 1 ' 56 ' 57 Kondo quenching of the 
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FIG. 3: Entropy S\ mp (T) vs T/ttF for systems in the MO 
regime. Plotted for common U/irT — 10 and t/nF = 5, with 
t'/7rT tuned to approach the transition [line (g), t' c /itT w 
2.45793] from the SC phase [t'/irT = 2.5, 2.459 and 2.45794 
for lines (a), (c) and (e)] and from the LM phase [t f /ttT — 2.4, 
2.457 and 2.45792 for lines (b), (d) and (f)]. The inset shows 
the excess charge, ni mp (T), for the same systems. 



TQD spin-^ is thus inoperative, with the TQD decou- 
pling from the leads on renormalization under reduction 
of the temperature/energy scale, leaving asymptotically 
a free local moment as T — » 0. The residual entropy at 
the corresponding LM fixed point is thus S^mp = ln(2), 
and leading irrelevant corrections to the fixed point are 
non-analytic. 56 ' 57 

Since the ground states for E A > (LM) and E A < 
(SC) cannot be continuously connected, a QPT is ex- 
pected around t' ~ U/A [see arrow (a) of Fig. 2]. This 
physical picture is confirmed directly in Fig. 3, where 
NRG results for the full TQD model are presented. 
Specifically, we show the TQD contribution to the en- 
tropy 5i mp (T) as a function of temperature, for sys- 
tems approaching the QPT between LM and SC phases. 
For T ^> l-E^J, the 3- and 4- electron states are quasi- 
degenerate, and a ln(3) entropy thus results. Flow to 
either the LM ground state with S'imp = ln(2) or the 
SC ground state with Si mp = occurs on the scale of 
Ta = \E' A \, which as such vanishes linearly with \t' — t' c \ 
as the transition at the critical t' c ~ Z7/4 is approached 
from either side (and which linearity is symptomatic of 
a parity-protected lev el- crossing QPT 29 ). 

The inset shows the excess charge 1 ' 5 due to the TQD, 
ni mp (T), for the same systems; confirming that the tran- 
sition at T = occurs between effective 3- and 4-electron 
TQD states. These results also show that precisely at the 
QPT (curve (g) in Fig. 3), the problem can be understood 
in terms of coexisting SC and LM fixed points - for ex- 
ample the T = charge n imp = 3| = ^t(3 x 2 + 4 x 1), 
is simply an average of the triple charge characteristic of 
the doubly-degenerate LM fixed point and the quadruple 
charge of the singly-degenerate SC fixed point. 
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C. Spin regime: strong interimpurity coupling 

More subtle behavior arises in the strongly-correlated 
case U ^> £, t',T, where each dot is essentially singly- 
occupied for temperatures T CU. However, the tunnel- 
couplings do generate effective exchange couplings be- 
tween the dots. This can be seen by projecting (SWT) 
the full Htqd onto the singly-occupied manifold of TQD 
states, perturbatively eliminating virtual excitations to 
2- and 4-electron TQD states to second-order in the 
tunnel-couplings t, t' and V. The resulting effective spin 
model is then given by 

H sp in = Pleads + J&1 ' (S 2 + S3) + J'S 2 • S3 + J K S 1 • s(O) 

(10) 

where S t is a spin-^ operator for dot i, s(0) is the con- 
duction electron spin density at the TQD as before, and 
the exchange couplings are given by 
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pJk = - Ft - (11) 



Two parity-distinct doublet ground states are then ob- 
tained for the isolated TQD, depending on whether t' > t 
(even parity) or t' < t (odd parity). In Ref. 29 we derived 
analytically an effective low-energy model to describe this 
lead-coupled TQD system for t 2 /U > T K but 
corresponding to regime (b) in Fig. 2. In either case, 
t' ^> t or t' <C £, the effective low-energy model is again 
a Kondo model Eq. 9, describing the residual coupling of 
the lowest TQD doublet state to the leads. Importantly 
however, j£ = —\Jk for t' <C t and </£ = +Jk for 
t' ^> t. The ferromagnetic Kondo effect 1 thus again arises 
for t' < t' c ~ t, with the lowest TQD doublet decoupling 
asymptotically at low energies. However, for t' > t f c the 
regular spin-^ antiferromagnetic Kondo effect 1 drives the 
system to a strong coupling state below the Kondo scale 
Tk ~ D exp(— 1/pJx), in which all TQD degrees of free- 
dom are quenched by the lead conduction electrons. The 
entropy is thus Si mp (0) = 0, corresponding to the Kondo 
singlet ground state. 1 ' 4 ' 5 On tuning t' through t' c ~ t, a 
level-crossing QPT occurs between LM and SC phases. 29 

One can in fact derive more generally an effective low- 
energy Kondo model of form Eq. 9, working directly with 
Htqd (i-e. without first projecting into the spin sec- 
tor, as above). The effective Kondo model is obtained 
simply by projecting the full Hamiltonian onto the re- 
duced TQD Hilbert space spanned by the lowest eigen- 
state of the numerically-diagonalized Htqd- For any 
t' < t' c ~ t <C U the odd-parity TQD doublet couples 
to the leads, while the even parity doublet couples for 
t f c < t f <C U. The QPT can of course also be realized 
by tuning t, keeping t' fixed (Fig. 2), and a numerical 
calculation of the effective coupling so obtained, j£, is 
shown in Fig. 4. The asymptotic values 29 J^/Jr = — | 
and +1 in the LM and SC phases are recovered in the 
singly-occupied limit; and the discontinuous change in 
the sign of is seen clearly at the point t' = £, due 
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FIG. 4: Effective coupling J^/Jk vs t/nT for the low-energy 
Kondo model Eq. 9 in the regime t,t ; <C U. Here t'/irT = 
10" 3 and U/ttT 8. 



to the level crossing of TQD states. NRG calculations 
on the full model (without any low-energy projections) 
indeed confirm that the transition occurs very close to 
t' = t. 



D. Spin regime: weak interimpurity coupling 

As suggested by arrow (c) in the phase diagram Fig. 2, 
a QPT between LM and SC phases arises even when the 
interdot couplings t and t' are very small. In this case 
however, the mechanism for local moment formation and 
Kondo screening is rather subtle. 

Since t, f, T <C U, the dots are still essentially singly- 
occupied, so the effective spin Hamiltonian Eq. 10 re- 
mains valid. In the limit where J = (arising when 
t = 0), dot 1 decouples from dots 2 and 3, and under- 
goes the regular Kondo effect with the lead to which it 
is coupled, being screened below T ~ Tk (with Tk the 
Kondo scale Tk ~ D exp(— 1/pJx) as above). 1 This must 
also in fact remain the case for small but non-vanishing 
J <C Tk, since the quenched Kondo singlet is already 
formed on the scale of Tk- However, this fixed point 
contains residual unquenched degrees of freedom corre- 
sponding to the spins on dots 2 and 3. In order to analyze 
the stability of this fixed point, one must thus determine 
whether there is an effective coupling of these degrees of 
freedom to the remaining Fermi liquid bath states of the 
lead. 

The effective exchange coupling acting directly be- 
tween dots 2 and 3 is given simply by J' = At' 2 /U 
(see Eqs. 10, 11), such that the local S = singlet 
state is lower in energy than the local S = 1 triplet 
by = E T — Es = J' • However, there is an addi- 
tional effective coupling between dots 2 and 3 due to an 
RKKY-type interaction mediated by the Kondo singlet 
formed between dot 1 and the leads. Virtual polariza- 
tion of this Kondo singlet is readily shown to generate 
an effective ferromagnetic contribution to the coupling 
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between the spins of dots 2 and 3; with second-order per- 
turbation theory in the coupling J within a Wilson chain 
formalism 21 ' 58 thereby yielding a renormalized singlet- 
triplet splitting, E' A = J' - \J 2 /T K (with A = (9(1) an 
undetermined positive constant). When E A > 0, the 2-3 
singlet thus lies lowest and decouples from the rest of the 
system. In consequence, at temperatures T <C \E A \ the 
entire system is in a singlet state, characteristic of the 
SC phase of regime (c) in Fig. 2. 

By contrast, when E A < the triplet formed between 
dots 2 and 3 lies lowest in energy. Thus in this case there 
are still residual TQD degrees of freedom at tempera- 
tures T <C \E A \. To determine the stability of this state, 
we must again consider the effective coupling between 
this triplet and the rest of the system. Within the Wil- 
son chain picture, it can be shown that the TQD triplet 
experiences an effective coupling to the remaining Fermi 
liquid bath states of the lead, mediated via the Kondo 
singlet formed with dot 1. The mechanism here is in fact 
completely analogous to that occurring in odd quantum 
dot chains with weak interdot coupling; 58 or in the asym- 
metric two-impurity Kondo problem. 21 For E A < 0, the 
effective low-energy model follows as 
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FIG. 5: Thermodynamics for the case of small interdot 
tunnel-couplings, t 2 /U <C Tk. The main panel shows the 
TQD contribution to entropy 5i mp (T) vs T/71T, while the in- 
set shows Tximp(T). Solid line for a system in the LM phase 
(E'a < 0) ; dashed line in the SC phase (E A > 0). U/ttT = 7, 
t/irT = 6 x 10" 3 and t' /ttT = 1 x 10" 3 (LM) or 6 x 10" 3 (SC). 



Hf =1 = H leads + J^S-s(0), 



(12) 



where S is now a spin-1 operator for the residual TQD 
state, i^ieads is the free conduction electron Hamiltonian 
with the Wilson chain 'zero orbital' removed due to the 
first-stage Kondo screening involving dot 1; and s(0) is 
its spin density at the TQD. The key result is that the 
effective coupling pj^ ~ t 2 /T K > is now antiferromag- 
netic, and the effective low-energy model Eq. 12 is thus 
a realization of the famous single-channel spin-1 Kondo 
model of Nozieres and Blandin 59 . 

The rich physical behavior of the spin-1 Kondo 
model 29 ' 56 ' 57 ' 59 61 is thus expected on the lowest energy 
scales. The antiferromagnetic coupling, J£ is renormal- 
ized upward on reduction of the temperature scale, re- 
sulting in quenching (or 'underscreening') of the S = 1 
'impurity' below T K - Dexp(-l/pJ*) (< <C 
Tk) 59 ; such that the ground state comprises a residual 
free spin-^ local moment, with leading RG-irrelevant cor- 
rections to the fixed point that are ferromagnetic (and 
non- an aly t ic ) . 5 6 ' 5 7 

In consequence, there is again a QPT on tuning f , 
separating a Kondo screened SC phase from a LM phase 
(regime (c) of Fig. 2). The transition itself is expected to 
occur at E A = J' — A J 2 /T K ~ 0, implying J' ~ J 2 /T K 
[or equivalently t' ~ t 2 /y/UTjc, see Fig. 2, arrow (c)]. 

To illustrate this complex behavior, Fig. 5 shows full 
NRG results for the TQD contribution to entropy and 
magnetic susceptibility for systems close to the transi- 
tion. Both systems have Tk/ttT & 10 -4 and similar 
\E a \/ttT & 10 -7 . The solid line corresponds to the LM 
phase, and the dashed line to the SC phase. On the tem- 
perature scale T ~ U, the dots become singly occupied, 

and so 5-^ = 3 ln(2) and Tximp = | corresponding to Solid lines are for E A > in the SC phase, approach- 
three free spins- 1 ,. As the temperature is reduced below ing progressively closely (from (d)-(a)) the transition at 



Tk, dot 1 is screened by the Kondo effect, leaving two 
quasi-degenerate spins on dots 2 and 3 (yielding thereby 
S imp = 21n(2) and T X i mp = f). For E' A > (dashed 
line), the residual TQD singlet state is lowest, and so 
Simp = and Tximp = when T <C \E A \. However, for 
E A < (solid line), the crossover is first to a free resid- 
ual TQD triplet state on the scale of \E A \; then below 
T ~ T K to the ultimate stable LM fixed point describing 
the underscreened S = 1 Kondo state, with Si mp = ln(2) 
and Tximp = \- 

Similar behavior is observed in dynamical quantities, 
such as the T = spectral function for dot 1, the en- 
ergy/frequency dependence of which is shown in Fig. 6. 
In the SC phase, the classic three-peak structure is ob- 
served. At high energies uj = ±U/2, the only spectral 
features are the Hubbard satellites, whose origin is sim- 
ply dot charge fluctuations. 1 At low energies \u\ ~ Tk, 
the narrow Kondo resonance is observed (see in par- 
ticular inset (A), dashed line); with the unitarity limit 
7rrDi(j(uj) = 1 reached at the Fermi level, uj = 0, being 
characteristic of the SC fixed point. In the LM phase 
however (solid lines), first-stage Kondo screening of dot 
1 on the scale of Tk is followed by second stage un- 
derscreening of the residual TQD triplet state, which 
as above is mediated via dot 1. This results in a slow 
crossover on the scale of T K to the LM fixed point, char- 
acterized by TiYDi a {uj) rsj 1/ In 2 \uj/T k \ behavior, such 
that TrTDicr^) = only at uj = 0, characteristic of the 
singular Fermi liquid, and reflecting the marginally irrel- 
evant corrections to the fixed point. 29 ' 56 ' 57 ' 60 ' 61 

Finally, in Fig. 7 we show the zero-bias conductance 
through the device as a full function of temperature. 
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FIG. 6: Single-particle spectrum of dot 1, nYDi a {uj) vs 
oj/ttF at T = 0, for systems with the same parameters as 
5. Insets (A) and (B) show detail of the low-energy be- 
havior on a logarithmic energy scale. The approach to the 
LM fixed point is characterized by marginally irrelevant loga- 
rithmic corrections, as highlighted in panel (B). However, for 
Tr <C \oj\ <C Tk the spectrum approaches the unitarity limit 
-KTDi a (<jo) — 1 due to the first-stage Kondo screening of dot 
1 [see panel (A)]. 



FIG. 7: Zero-bias conductance through the TQD, G c (T)/Gq 
vs T/T£ jC , for systems closely approaching the transition. 
T Kc is the crossover scale upon which the transition fixed 
point is reached at t = t c . Plotted for U/ttT = 4 and 
t'/irT = 10" 2 , varying t/vrr = t c ±\T K , with A = 10" 7 , 10" 5 , 
10 -3 and 10 _1 for lines (a)-(d) respectively. Solid lines for 
systems in the SC phase; dashed lines for the LM phase. Tk 
is the common scale for the first-stage Kondo effect involving 
dot 1, as indicated on the plot. 



E f A = 0; while dashed lines are for E' A < in the LM 
phase. All systems have a common scale Tr-, and the 
data are rescaled in terms of T/Tj Cc (with T£ c chosen 
arbitrarily as the T for which G c (T)/Go = |). The 
S = 7r/2 phase shift in the Kondo screened SC phase 1 
implies a unitarity T = conductance G c (0)/Go = 1 
from Eq. 7; while S = in the LM phase 60 gives van- 
ishing conductance, 60 G c (0)/Go = 0. The behavior at 
the transition fixed point itself can again be understood 
here in terms of coexisting LM and SC states; yielding a 
T = conductance at the transition of G c (0)/Go = |. 
The full temperature dependence is of course rich, reflect- 
ing with decreasing T renormalization group flow to the 
Kondo SC fixed point on the scale of T ~ Tr-, followed 
by flow to the transition fixed point at T^ c ; and finally 
flow to the stable LM or SC fixed points describing the 
true ground state, on the scale T ~ \E' A \. 



IV. RESULTS: DISTORTIONS 

We now turn to the more general case of broken parity 
symmetry, as occurs due to distortions of the triangu- 
lar TQD structure (Fig. 1). The presence of symmetry- 
breaking perturbations, such as ti2 ^ £13, of course pre- 
cludes labeling states by a parity quantum number, since 
[H, P23] 7^ 0. In the symmetric case considered above, 
the QPT between LM and SC phases was characterized 
by a level crossing between parity-distinct states. We 
show below that the parity-broken model still supports 
LM and SC phases, but that the level crossing transition 



becomes instead a QPT of Kosterlitz-Thouless form. 49 
The schematic phase diagram, Fig. 2, still applies in the 
general case, but with the phase boundary now to be 
understood as a line of SC critical end points, with the 
Kondo scale vanishing as the transition is approached 
from the SC phase. 



A. Effective models 

We focus here on the strongly correlated case of pri- 
mary interest, U ^> T. At temperatures T <C U, each 
dot becomes essentially singly-occupied, and an effective 
TQD spin model analogous to Eq. 10 can similarly be 
derived via a SWT upon perturbative elimination of vir- 
tual excitations to 2- and 4-electron TQD states. For ease 
of comparison with the parity-symmetric case, we retain 
tis = t and £23 = t', defining the asymmetry in terms 
of the parameter x = £12 — £13. The resultant effective 
low-energy model then has the form 

#spin = ^leads + ^S-L • (S 2 + S 3 ) + J'S 2 • S3 + J K S 1 • s(0)+^S 1 • S 2 

(13) 

where the effective exchange couplings J, J' and Jk are 
given in Eq. 11; and where the asymmetry enters only 
through the final term, with 

5 = 4^±^1. (14) 

The low-energy manifold of the isolated TQD again 
comprises a pair of doublet states, denoted \A\S Z ) and 
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FIG. 8: Avoided crossing of the lowest TQD doublet states in 
the asymmetric case shown in the inset, with 5/ J — 0, 0.2, 0.4 
and 0.6 for lines (a)-(d) respectively. Main panel shows the 
effective coupling of the lowest doublet to the leads Jk^e/ Jk 
vs J'/J' c , for U/ttT = 10, t/irT = 0.1 and 6/irT = KT 4 , 
obtained from Eq. 19 (solid line); and extracted from NRG 
data in the SC phase via Tk ~ D exp(— l/pJK,es) (points). 



\B;S Z ) (with an essentially irrelevant spin quartet oc- 
curring at higher energies). For any degree of asymme- 
try, these states may be expressed in terms of even- and 
odd-parity doublets |±; S z ) (themselves obtained in the 
symmetric limit x = = 5 considered above). They are 
given by 



\A;S*)= + a\-,S*)+l3\+;S z ) 
\B;S z }=-f3\-,S z )+a\+;S z } 



where 



1-/3 2 



1 + 



J-J' + \5 



(15a) 
(15b) 



• (16) 



The energy difference between these isolated TQD dou- 
blet states is, 



E A = E TQD (B) - E TQD (A) = itf)2+ 3 5 2 

(17) 

such that E/\ > for all finite values of J, J' and 5 (or 
t, t' and x). The effect of asymmetry is thus to turn 
the level crossing of TQD states at 5 = to an avoided 
crossing when 5^0. There is in other words a single, 
unique doublet TQD ground state (the \A;S Z ) state) for 
all model parameters in the asymmetric case. This is 
illustrated in the inset to Fig. 8, where we show the en- 
ergies of the TQD doublet states as a function of J' for 
different values of the asymmetry 5 (the solid line (a) is 
for the symmetric case 5 = 0). 

The implication of Eq. 17 is that below the tem- 
perature/energy scale only the TQD doublet state 



| A; S z ) is accessible. In consequence one can derive a low- 
energy effective model, valid for T <C £?a, by projecting 
the spin model Eq. 13 onto the ground state TQD dou- 
blet manifold spanned by \A;S Z ). This can be done to 
first order in the Kondo coupling J^, leading to 



H eff =H leads + J K }2 \A ] S z )(A ] S z \S 1 \A ] S z ')(A ] S z '\.s(0) 

S z ,S zf 

=H\ ea ds + JK,effS ' S(0) , 

(18) 

where S is a spin-^ operator for the doublet TQD state 
\A\ S z ); and where the effective exchange coupling to the 
leads follows as 



JjK,eff= Jk 1 



-a 



(19) 



From Eq. 16, we immediately find that Jx,eff > is anti- 
ferromagnetic for J' > J, but is ferromagnetic, Jx,eff < 
0, for J' < J (and one also recovers asymptotically the 
results of the parity-symmetric case, JK,eff — >* —\Jr, for 
J > J', 5; and J K ,eff -> +Jk for J' > J,£). The effec- 
tive coupling as a function of J' / J' c is plotted in Fig. 8 as 
the solid line; points are JK,eff extracted from full NRG 
results, and show very good agreement. 

Importantly, the antiferromagnetic effective Kondo 
coupling vanishes continuously as (J' — J) — » + [or 
(t' - t) 0+] for any \5\ > 0, 



J K,eff-- 



Jk 
25 



(J' - J) + G(J' - Jf 



(20) 



(in contrast to the discontinuous change in the effective 
Kondo coupling associated with a level crossing QPT, see 
e.g. Fig. 4). In consequence, in the full model we expect 
a Kondo SC phase for t' > t' c {p± t), with TQD degrees of 
freedom entirely quenched 1 below a Kondo scale Tk ~ 
Z}exp(— 1/pJk eff), which itself vanishes continuously as 
(j> _ j*) _> 0+' [or (tf-^) -►()+)]: 



Tk ~-Dexp 
~£> exp 



25 \ 1 
Jhc) J'-J' c 
irUx(2t + x)\ 1 



f - 1' 



(21) 



The form of the vanishing Kondo scale as the transition 
is approached, Eq. 21, is that arising for a Kosterlitz- 
Thouless 49 (KT) transition. The point t' = t' c can thus 
be understood as the critical end point of a line of Kondo 
SC fixed points. For t' < t' c , there is no low-energy scale; 
and the residual TQD doublet degree of freedom remains 
free down to T — >> at this LM fixed point. 



B. Kosterlitz-Thouless transition 

This is explored further in Figs. 9 and 10, containing 
NRG results for the full TQD model as the KT transition 
is approached from the SC phase, t r > t' c . 





FIG. 9: Thermodynamics for the parity-broken model. TQD 
contribution to entropy S\ mp (T) (main panel) and magnetic 
susceptibility Tx imp (T) (inset) vs T/tiT. Plotted for U/ttT = 
10, t/vrr = 0.1 and S/ttT = 10~ 4 , varing t'/irT = 9.8 x 10~ 2 , 
9.78 x 10" 2 , 9.76 x 10" 2 , 9.74 x 10" 2 and 9.72 x 10" 2 to 
approach progressively the transition. Kondo scales Tk are 
indicated by vertical arrows, and vanish continuously as the 
transition is approached. 



Fig. 9 shows thermodynamic quantities Si mp and 
^Ximp as a function of temperature. On the scale T ~ U, 
the dots become singly-occupied, and so the entropy 
drops to ^imp = 31n(2), with a corresponding tripled 
Curie law susceptibility, Tximp = f • The quartet TQD 
state becomes inaccessible on the scale T ~ J leaving 
the pair of quasi-degenerate TQD doublet states (and 
hence S[ mp = In (4)). The higher lying doublet is in turn 
projected out on the scale of This is the LM fixed 
point, with S- imv = m(2) and Tximp = \- The lowest 
TQD doublet is then screened by the Kondo effect on 
the scale of Tk, which itself vanishes as the transition 
is approached (as evident directly from the figure). At 
the SC fixed point, TQD degrees of freedom are entirely 
quenched, giving S imp = and Tx imp = 0. 

Similar behavior is observed in the T = spectral func- 
tion of dot 1 shown in Fig. 10, together (inset) with the 
T-dependence of the zero-bias conductance through the 
TQD. On the energy scale \uj\ ~ E&, the higher energy 
TQD doublet state \B;S Z ) is projected out, terminating 
the renormalization of its coupling to the leads. Instead, 
the lowest TQD \A; S z ) doublet flows to strong cou- 
pling, with the characteristic narrow Kondo resonance 
evident in the spectrum on the scale of \uj\ ~ Tk\ and 
reaching the unitarity limit i\YD\ a {uj = 0) = 1 that is 
characteristic 1 of the strong coupling fixed point. The 
Kondo resonance narrows progressively as Tk dimin- 
ishes and the QPT is approached from the SC phase; 
vanishing continuously 'on the spot' at the transition 
itself, such that in (and throughout) the LM phase 
itYDiv^uo = 0) = 0. Corresponding behavior is nat- 
urally seen in the T-dependence of the conductance, 



FIG. 10: Single-particle spectrum of dot 1, nYD\ a (uj) vs oj/ttT 
(on a log-scale) at T = 0, for systems with the same parame- 
ters as 9. Vertical arrows indicate the Kondo scales Tk- Inset 
shows the temperature-dependence of the zero-bias conduc- 
tance G c (T)/Go for the same systems. 



which is determined by the local spectrum as in Eq. (6). 
Throughout the SC phase, it is likewise enhanced at low- 
temperatures T < Tk clue to the Kondo effect, with 
G C (T = 0)/G = 1 in all cases. In the LM phase by 
contrast, where 7iTDi a {uj = 0) = 0, the T = conduc- 
tance vanishes - the conductance thus changing discon- 
tinuously as the transition is crossed, which appears to be 
a rather general signature of KT transitions in quantum 
dot and related systems. 60 ' 62 65 



V. SUMMARY AND DISCUSSION 

A compact cluster of three quantum impurities, hy- 
bridizing apically with a single channel of host conduction 
electrons, has been shown to exhibit a rich range of phys- 
ical behavior. Such a model may describe apex-coupled 
impurity trimers on metallic surfaces, or semiconductor 
TQD devices. 

In the parity-symmetric case where two of the impuri- 
ties/dots are in equivalent local environments, a level- 
crossing quantum phase transition separates a Kondo 
screened strong coupling (SC) phase and a free local mo- 
ment (LM) phase. For strong interdot tunnel-coupling 
(and hence weak interactions), the transition in essence 
occurs between states with either 3 or 4 electrons occu- 
pying TQD molecular orbit als. For weaker couplings and 
strong dot electron correlations, sites of the TQD are es- 
sentially singly-occupied. Within this TQD spin-regime, 
the QPT between SC and LM phases again arises; al- 
though different mechanisms are uncovered, depending 
on the relative interdot tunneling strengths. For weak 
interdot couplings in particular, the Kondo effect screens 
a single dot, with the residual TQD degrees of freedom 
forming either a local singlet (SC) or an underscreened 
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spin-1 Kondo state (LM). When parity symmetry is bro- 
ken by distortions of the triangular TQD structure, we 
show that LM and SC phases are still supported — with 
the transition between them now of Kosterlitz-Thouless 
type, such that the low-energy Kondo scale vanishes con- 
tinuously as the transition is approached from the SC 
phase. 

Given the ubiquity of SC and LM phases, an obvious 
question arises: what characterizes - and distinguishes - 
such phases more generally? One answer is to note that 
the T = zero-bias conductance is pinned in the singly- 
occupied case to a unitarity value Gf c (T = 0)/Go = 1 

in the Kondo SC phase, but vanishes in the LM phase, 
G LM( T = Q y GQ = The Qrigin of thijg can be tmced 

via Eq. 7 to the scattering phase shift, 5, experienced by 
conduction electrons in either phase. 

By a straightforward extension to the TQD model of 
the analysis given in Ref.60, a Friedel-Luttinger sum 
rule 60 ' 66 ' 67 can be shown to relate exactly the phase shift 
to the excess charge (ni mp (T = 0)) via S = fni mp (T = 
0) + II; where the Luttinger integral II is given by 

I L = ImTr f du^^G^u) (22) 



(with G a (uj) and T> a (uj) the 3x3 matrices for the lead- 
coupled TQD Green functions and self-energies respec- 
tively). We find that the SC and LM phases are each 
characterized by a distinct value for the 'topological' 
quantity II' from extensive NRG calculations of the 
Luttinger integral, and regardless of the bare underlying 
model parameters, we find II = to be characteristic 
of the SC phase (as indeed is well known for any Fermi 
liquid phase 66 ), and \Il\ = § to be equally characteristic 
of the LM phase. The latter is also precisely as found 
throughout the entire LM phase of a correlated 2-level 
quantum dot model 60 , suggesting that the Luttinger in- 
tegral is indeed a universal characteristic of a LM phase 
- albeit that the fundamental reasons for this await an 
answer. 



Acknowledgments 

We are grateful for funding from EPSRC (UK) through 
grant EP/1032487/1. 



1 A. C. Hewson, The Kondo Problem to Heavy Fermions 
(Cambridge University Press, Cambridge, 1993). 

2 J. Kondo, Prog. Theor. Phys. 32, 37 (1964). 

3 P. W. Anderson, Phys. Rev. 124, 41 (1961). 

4 K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

5 H. R. Krishnamurthy, J. W. Wilkins, and K. G. Wilson, 
Phys. Rev. B 21, 1003, 1044 (1980). 

6 D. Withoff and E. Fradkin, Phys. Rev. Lett. 64, 1835 
(1990). 

7 C. Gonzalez-Buxton and K. Ingersent, Phys. Rev. B 54, 
15614 (1996). 

8 D. E. Logan and M. T. Glossop, J. Phys.: Condens. Matter 
12, 985 (2000). 

9 R. Bulla, M. T. Glossop, D. E. Logan, and T. Pruschke, J. 
Phys.: Condens. Matter 12, 4899 (2000). 

10 M. T. Glossop and D. E. Logan, J. Phys.: Condens. Matter 
15, 7519 (2003). 

11 L. Fritz, and M. Vojta. arXiv: 1208.3113 (2012). 

12 A. K. Mitchell, and L. Fritz. arXiv: 1212 . 2631 (2012). 

13 L. Fritz and M. Vojta, Phys. Rev. B 72, 212510 (2005). 

14 M. V. A. K. Mitchell, D. Schurucht and L. Fritz, Phys. 
Rev. B 87, 075430 (2013). 

15 R. Zitko, Phys. Rev. B 81, 241414(R) (2010). 

16 B. A. Jones, C. M. Varma, and J. W. Wilkins, Phys. Rev. 
Lett. 61, 125 (1988); B. A. Jones, Physica B (Amsterdam) 
171, 53 (1991). 

17 I. Affleck and A. W. W. Ludwig, Phys. Rev. Lett. 68, 1046 
(1992); I. Affleck, A. W. W. Ludwig and B. A. Jones, Phys. 
Rev. B 52, 9528 (1995). 

18 F. W. Jayatilaka, M. R. Galpin, and D. E. Logan, Phys. 
Rev. B 84, 115111 (2011). 

19 E. Sela, A. K. Mitchell, and L. Fritz, Phys. Rev. Lett. 106, 
147202 (2011). 

20 A. K. Mitchell and E. Sela, Phys. Rev. B 85, 235127 



(2012). 

21 A. K. Mitchell, E. Sela, and D. E. Logan, Phys. Rev. Lett. 
108, 086405 (2012). 

22 C. Jayaprakash, H. R. Krishnamurthy, and J. W. Wilkins, 
Phys. Rev. Lett. 47, 737 (1981). 

23 W. G. van der Wiel et al., Rev. Mod. Phys. 75, 1 (2003). 

24 A. M. Chang and J. C. Chen, Rep. Prog. Phys. 72, 096501 

(2009) . 

25 J. Merino, L. Borda, and P. Simon, Europhys. Lett. 85, 
47002 (2009). 

26 S. J. Chorley, M. R. Galpin, F. W. Jayatilaka, C. G. Smith, 
D. E. Logan, and M. R. Buitelaar, Phys. Rev. Lett. 109, 
156804 (2012). 

27 T. Jamneala, V. Madhavan, and M. F. Crommie, Phys. 
Rev. Lett. 87, 256804 (2001). 

28 K. Ingersent, A. W. W. Ludwig, and I. Affleck, Phys. Rev. 
Lett. 95, 257204 (2005). 

29 A. K. Mitchell, T. F. Jarrold, and D. E. Logan, Phys. Rev. 
B 79, 085124 (2009). 

30 A. K. Mitchell and D. E. Logan, Phys. Rev. B 81, 075126 

(2010) . 

31 T. Uchihashi, J. Zhang, J. Kroger, and R. Berndt, Phys. 
Rev. B 78, 033402 (2008). 

32 B. Lazarovits, P. Simon, G. Zarand, and L. Szunyogh, 
Phys. Rev. Lett. 95, 077202 (2005). 

33 W. Z. Wang, Phys. Rev. B 76, 115114 (2007). 

34 Y. C. Xiong, J. Huang, and W.-Z. Wang, J. Phys.: Con- 
dens. Matter 24, 455604 (2012). 

35 R. Zitko and J. Bonca, Phys. Rev. B 77, 245112 (2008). 

36 A. Oguri, S. Amaha, Y. Nishikawa, T. Numata, M. Shi- 
mamoto, A. C. Hewson, and S. Tarucha, Phys. Rev. B 83, 
205304 (2011). 

37 T. Numata, T. Nisikawa, A. Oguri, and A. C. Hewson, 
Phys. Rev. B 80, 155330 (2009). 



12 



38 Y. B. Kudasov and V. M. Uzdin, Phys. Rev. Lett. 89, 
276802 (2002). 

39 H. J. Gotsis, N. Kioussis, and D. A. Papaconstantopoulos, 
Phys. Rev. B 73, 014436 (2006). 

40 A. Antal, B. Lazarovits, L. Szunyogh, B. Ujfalussy, and 
P. Weinberger, Phys. Rev. B 77, 174429 (2008). 

41 D. Schroer, A. D. Greentree, L. Gaudreau, K. Eberl, 
L. C. L. Hollenberg, J. P. Kotthaus, and S. Ludwig, Phys. 
Rev. B 76, 075306 (2007). 

42 A. Vidan, R. Westervelt, M. Stopa, M. Hanson, and 
A. Gossard, J. Supercond. Incorp. Novel Magn. 18, 223 
(2005). 

43 M. C. Rogge and R. J. Haug, Phys. Rev. B 77, 193306 

(2008) . 

44 L. Gaudreau, A. Kam, G. Granger, S. A. Studenikin, P. Za- 
wadzki, and A. S. Sachrajda, Appl. Phys. Lett. 95, 193101 

(2009) . 

45 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. 
Abusch-Magder, U. Meirav, and M. A. Kastner, Nature 
391, 156 (1998); S. M. Cronenwett, T. H. Oosterkamp, 
and L. P. Kouwenhoven, Science 540, 281 (1998). 

46 M. A. Kastner, Physics Today 46, 24 (1993). 

47 H. Jeong, A. M. Chang, and M. R. Melloch, Science 293, 
2221 (2001). 

48 R. H. Blick, R. J. Haug, J. Weis, D. Pfannkuche, K. v. 
Klitzing, and K. Eberl, Phys. Rev. B 53, 7899 (1996). 

49 J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 
(1973). 

50 R. Bulla, T. Costi, and T. Pruschke, Rev. Mod. Phys. 80, 
395 (2008). 

51 R. Peters, T. Pruschke, and F. B. Anders, Phys. Rev. B 
74, 245114 (2006). 



52 A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 
076402 (2007). 

53 F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 
(2005). 

54 Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992). 

55 J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 
(1966). 

56 P. Mehta, N. Andrei, P. Coleman, L. Borda, and 
G. Zarand, Phys. Rev. B 72, 014430 (2005). 

57 W. Roller, A. C. Hewson, and D. Meyer, Phys. Rev. B 72, 
045117 (2005). 

58 A. K. Mitchell, D. E. Logan, and H. R. Krishnamurthy, 
Phys. Rev. B 84, 035119 (2011). 

59 P. Nozieres and A. Blandin, J. Phys. (Paris) 41, 193 (1980). 

60 D. E. Logan, C. J. Wright, and M. R. Galpin, Phys. Rev. 
B 80, 125117 (2009). 

61 C. J. Wright, M. R. Galpin, and D. E. Logan, Phys. Rev. 
B 84, 115308 (2011). 

62 W. Hofstetter and H. Schoeller, Phys. Rev. Lett. 88, 
016803 (2002). 

63 M. Garst, S. Kehrein, T. Pruschke, A. Rosch, and M. Vo- 
jta, Phys. Rev. B 69, 214413 (2004). 

64 M. R. Galpin, D. E. Logan, and H. R. Krishnamurthy, 
Phys. Rev. Lett. 94, 186406 (2005). 

65 M. R. Galpin, D. E. Logan, and H. R. Krishnamurthy, J. 
Phys.: Condens. Matter 18, 6545 (2006). 

66 J. M. Luttinger and J. C. Ward, Phys. Rev. 118, 1417 
(1960). 

67 D. C. Langreth, Phys. Rev. 150, 1516 (1966). 



